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ABSTRACT 

Context. This methodological paper is part of a short series dedicated to the long-standing astronomical problem of 
de-projecting the bi-dimensional, apparent morphology of a three-dimensional distribution of gas. 

Aims. We focus on the quantification and spatial recovery of turbulent motions in planetary nebulae (and other classes 
of expanding nebulae) by means of long-slit echellograms over a wide spectral range. 

Methods. We introduce some basic theoretical notions, discuss the observational methodology, and develop an accurate 
procedure disentangling all broadening components (instrumental resolution, thermal motions, turbulence, gradient of 
the expansion velocity, and fine structure of hydrogen-like ions) of the velocity profile in all spatial positions of each 
spectral image. This allows us to extract random, non-thermal motions at unprecedented accuracy, and to map them 
in 1-, 2- and 3-dimensions. 

Results. We discuss general and specific applications of the method. We present the solution to practical problems in 
the multi-dimensional turbulence-analysis of a testing-planetary nebula (NGC 7009), using the three-step procedure 
(spatio- kinematics, tomography, and 3-D rendering) developed at the Astronomical Observatory of Padua (Italy). In 
addition, we introduce an observational paradigm valid for all spectroscopic parameters in all classes of expanding 
nebulae. 

Conclusions. Unsteady, chaotic motions at a local scale constitute a fundamental (although elusive) kinematical pa- 
rameter of each planetary nebula, providing deep insights on its different shaping agents and mechanisms, and on 
their mutual interaction. The detailed study of turbulence, its stratification within a target and (possible) systematic 
variation among different sub-classes of planetary nebulae deserve long-slit, multi-position angle, wide-spectral range 
echellograms containing emissions at low-, medium-, and high-ionization, to be analyzed pixel-to-pixel with a straight- 
forward and versatile methodology, extracting all the physical information (fiux, kinematics, electron temperature and 
density, ionic and chemical abundances, etc.) stored in each frame at best. 

Key words, planetary nebulae: general- ISM: kinematics and dynamics- ISM: structure 
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els (jPerinotto et al. I 2004b; ISchonberner et al. I 2005) still 
partly fail in reproducing the positive expansion velocity 
gradient of the main shell (the Wilson's law; Wilson 1 1950; 
IWeedman~l 1968), and the gas distribution observed in a 
representative sam ple of carefully de-projected "true" PNe 
(|Sabbadin et al. 1 2006 and references therein). Such model- 
PN vs. true-PN discrepancies suggest that the former tend 
to systematically overestimate the relative importance of 
wind interaction (as supported by a recent revision down- 
ward of mass-loss rates in PN central sta rs, du e to clumping 
of the fast stellar wind: iKudritzki et al. l 2006; |Prinja et al. | 
2007). 

An even deeper uncertainty surrounds the turbu- 
lence of expanding g fundamental parameter mea- 
suring local, random, unsteady motions due to non-linear 
energy transfer. On the one hand, 2-D hydrodynami- 
cal simulations of wind interaction (Frank & Mellcma 
1994: IVishiiiacI 1994: fMellema . 1997: ,Dwarkadas fc Balick , 
1998; IGarcia-Segura et al. | 1999, 2006) predict strati- 
fied turbulent motions peaking within shocked regions 
and/or ionization fronts characterized by the growth of 
Kelvin-Hclmholtz, Vishniac or Ray leigh- Taylor instabili- 
ties. On the other hand, practical determinations based on 
high-r esolution spect ra ([Guerrero et al. 1 1998; I Acker et al. I 
2002; [Gesicki et ap 2003 , 2006; [Gesicki fc Zijlstra | 2003, 
2007; and Medi na et aTl 2006) provide a unique, average 
turbulence- value for each nebula, ranging from km s~^ to 
20 km s-^. 

We aim at quantifying and mapping turbulent motions 
in PNe by means of long-slit, wide spectral range, very- 
high spectral and spatial resolution echellograms at several 
position angles (PA), reduced and analyzed according to 
spatio-kinematical, tomographic and 3-D recovery meth- 
ods developed at the Astronomical Observatory of Padua 
(jSabbadin et al. 1 2000. 2006, and references therein). 

This introductory, methodological paper (a) dissects all 
broadening agents of a velocity profile; (b) describes an 
original procedure mapping turbulent motions in 1-, 2-, and 
3-dimensions; (c) illustrates the solution to practical prob- 
lems of turbulence-determination in a testing-PN; and (d) 
introduces an observational paradigm valid for all types of 
expanding nebulae. 

In a forthcoming paper, we (i) will present multi- 
dimensional turbulence maps for a representative sample of 
PNe in different evolutionary phases; and (ii) will compare 
these observational results with expectations from hydro- 
dynamical simulations and theoretical models. 

2. Basic considerations 

A PN consists of (large- and small-scale) interconnected 
sub-systems with different morphology, physical conditions, 
and kinematics. Each elementary volume within the radial 
slice of nebula selected by a spectrograph slit is character- 
ized by 

1) a regular expansion, and 

2) random (i. e., thermal and turbulent) motions. 

Thanks to the large stratification of the radiation and the 
kinematics, each ionic species generates the characteristic 
bowed-shaped spectral image sketched in Fig. 1. 

Two useful and complementary reference lines of a 
spectral image ([Sabbadin et al. 1 2000. 2004, and references 
therein) are: 




Fig. 1. Nebular projection on the sky vs. high-dispersion 
spectral image. The radial slice of the Saturn Nebula in- 
tercepted by a spectrograph slit crossing the central star 
generates a series of high-resolution spectral images; each 
spectral image contains kinematical and spatial information 
on the parent ionic species. 

- the central-star-pixel-line (cspl), common to all PA, 

representing the emission of the matter projected at the 
apparent position of the central star, whose motion is 
purely radial; 

- the zero- velocity-pixel-column (zvpc), giving the spa- 

tial intensity profile of the tangentially moving gas at 
the systemic radial velocity. 

As illustrated in Fig. 2, the cspl of a spectral image 
contains the kinematical information: the intrinsic velocity 
profile of each nebular sub-system, convolved with a se- 
ries of broadening agents (instrumental resolution, thermal 
motions, turbulence, and fine-structure of hydrogen-like 
ions), produces an observed, "Gaussianized" velocity pro- 
file, whose peak position and fuU-width-at-half-maximum 
(FWHM) of the blue- and red-shifted components give the 
radial velocity and the velocity spread of approaching and 
receding layers, respectively. 

Analogously, the zvpc of a spectral image contains the 
spatial information: the intrinsic spatial profile of each neb- 
ular sub-system, convolved for seeing-|-guiding effects, pro- 
duces a "Gaussianized" spatial profile, whose peak separa- 
tion and FWHM give the angular distance and the thick- 
ness, respectively, of the corresponding ionic layers in the 
plane of the sky. 

Let us consider the cspl, and focus on the FWHM 
of blue- and red-shifted components (VFf ota/ (blue) and 
W^tota/(red), respectively). For a typical nebular sub-system, 
the "Gaussianization" of the cspl is very effective, and the 
value of Wtotai is the sum (in quadrature) of different broad- 
ening contributions, namely: 

(A) The gradient of the expansion velocity, corre- 
sponding to the range in radial velocity of each emit- 
ting layer, which is usually assimilable to a Gaussian 
profile with Wgrad- The value of Wgrad changes from 
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Fig. 2. Anatomy of a high-resolution spectral image, show- 
ing the location and flux distribution of both the central- 
star-pixel-line (cspl) and the zero-velocity-pixel-column 
(zvpc). More details are in the text. 



one ion to another, and can be estimated from the ex- 
pansion law and the spatio-kincmatical reconstruction. 
A common case concerns an ellipsoidal PN seen at ran- 
dom orientations and expanding with T4xp=Axr"; we 
have Wgrad=Ax Ar", where the angular thickness Ar" 
is given by the zvpc profile. Typical Wgrad values in PN 
main-shells are 10 to 15 km s^^ for H I, 7 to 10 km s^^ 
for He I and He II, and 4 to 6 km s~^ for ionic species 
of heavier elements. 

(B) Instrumental resolution, corresponding to a 
Gaussian profile with FWHM {Wmstr), is given by the 
spectral resolution (A/AA) of the spectrograph; a typi- 
cal value of A/AA=50 000 provides Wi„str=6.0 km s~^. 

(C) Thermal motions, generating a Gaussian dis- 
tribution with T/Ffher=(8xkxTcXln 2/m)°-^= 

km s~^, where Tp is the local 



0.214xT0-5xm 



-0.5 



electron temperature, and m is the atomic weight of 
the element. Thus, for Tc=10'' K, Wther amounts to 
21.4 km s~^ for hydrogen, 10.7 km s~^ for helium, 5.8 
km s^^ for nitrogen, and 5.4 km for oxygen. 
(D) Turbulence, i. e., non-thermal, chaotic motions at a 
local scale, is characterized by Wt„rb0 Due to the intrin- 
sically dissipative nature of turbulence (and in agree- 



ment with hydrodynamical simulations; e. g. JMellema I 
1997; IDwarkadas fc Balick 1 1998; |Garcia-Seg ura et al. | 
1999), Wturb varies across the nebula, peaking in 
shocked layers and ionization fronts. Thus, different 
ionic species are affected by different turbulent motions, 
and, on the analogy of the radiation and the kinemat- 
ics, we can introduce the concept of "stratification of 
the turbulence" . This means that a thorough investi- 
gation of the gas turbulence needs spectra containing 
emissions at high-, medium-, and low-ionization. 
(E) The fine-structure of hydrogen-like ions, whose 
levels with the same principal quantum number have 
small energy differences due to the spin-orbit cou- 
pling. We have performed a detailed analysis for 
Case B of IBaker fc Menzel I (1938), and a Z-levels 
population proportional to their statistical weights 
dStorey fc Hummer | 1995), obtaining M^/,=7.7(±0.3) 
km s~^ and 5.9(±0.3) km s^^ for the seven fine- 
structure components of H I Ha and H/3, respectively 
(also see |Clegg etlFI 1999); 10.1(±0.3) km s'^ for A 
4686 A of He II (thirteen fine-structure components); 
4.2(±0.3) km s^^ for A 6560 A of He II (nineteen fine- 
structure components); and 3.0(±0.3) km s~^ for A 5876 
A of He I (six fine-structure components). Moreover, the 
fine-structure profile of He II A 4686 A presents a weak 
- although showy - blue-shifted tail, and the one of He 
I A 5876 A a modest red-shifted tail. The values of Wjs 
are nearly independent on To for 5000 K<To<20000 K, 
and on the electron density {N^) for 10^ cm~^ <A'o<10^ 
cm 

For a typical nebular sub-system, the net turbulence- 
broadening in each ionic species is 



Wturb = [W^otal - 



grad 



-w- 



2 

■ her 



instr 



(1) 



The application of Eq. (1) to the cspl of emissions 
at low-, mean-, and high-ionization, combined with de- 
projection criteria introduced bv ISabbadin et al~l (2005), 
provides the overall 1-D turbulence-map of the nebula. 

To be noticed: beyond the cspl, Eq. (1) applies to each 
pixel-line along the dispersion, thus disentangling all broad- 
ening agents of the velocity profile at all "latitudes" of a 
spectral image. A different-latitude-pixel-line is identified 
by its angular distance (d") from the cspl; on the analogy 
of the latter, in the following it will be named "dlpl" . 

Thus, 

I - The application of Eq. (1) to all dlpl of a spectral im- 

age allows us to map the turbulence of the parent ion 
within the entire de-projected nebular slice covered by 
a spectrograph slit (ionic 2-D turbulence-map); 

II - The combination of all ionic 2-D turbulence-maps at a 
given PA provides the overall 2-D turbulence- map 
at the corresponding PA; and 

HI - By assembling overall 2-D turbulence-maps at all ob- 
served PA (for example, through the 3-D rendering pro- 
cedure introduced by Ragazzoni et al. 2001), we infer 
the overall 3-D turbulence-map of the entire neb- 
ula. 

Summing up, the determination of the local turbulence in 
1-, 2-, and 3-dimensions passes through the constraining of 



^ We "assume" a Gaussian probability distribution for the tur- 
bulence velocity, although some authors suggest a non-Gaussian 



behavior for vorticity, energy dissipation and passive scalars (for 
a review, see piiie green fc Scalo | 2004 ). 
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sinergistic broadening effects of both the instrumentation 
and nebular gas parameters (the latter include atomic prop- 
erties, physical conditions, kinematics and ionization struc- 
ture). This implies that an adequate instrumental spectral 
resolution >50000) is the "conditio sine qua non" 

for a quantitative study of turbulent motions in PNe, to be 
combined with a detailed spatio-kinematical reconstruction 
at high-, medium-, and low-ionization. 

In addition, large, extended nebular sub-systems (e. g., 
optically-thin outer shells) with a broad, flat-topped intrin- 
sic velocity profile departing from a Gaussian distribution, 
provide less reliable turbulence-values than localized, com- 
pact sub-systems with a sharp intrinsic velocity profile (e. 
g., the main shell, condensations and knots, etc.). Last, 
recombination lines of H I, He I and He H are less reli- 
able turbulence-diagnostics than coUisionally-excited lines 
of heavier elements. 

All this is illustrated in Fig. 3, showing [O HI] and H 
I spectral images of a random seen, regularly expanding 
(T4xpCxr) model-spheroid denser at the equator than at the 
poles, covered along the apparent major axis. We adopt 
three representative values of the instrumental resolution 
(A/AA=100 000, 50 000, and 25 000). Moreover, we assume 
To^lO"* K, Wturb—W km s^^, seeing=1.0", and instrumen- 
tal spatial and spectral scales of 0.2G"px~^ and 1.5 km s^^ 
px~^, respectively. Different acronyms used in the current 
discussion (cspl, dlpl, zvpc etc.) are indicated in Fig. 3, to- 
gether with straight lines (a and b) containing polar and 
equatorial axes of the nebular slice, and their projections 
(a' and b') on a spectral image. A ccording to the second de- 
projection rule introduced bv i Sabbadin et al71 (2005), x-y 
coordinates of the nebular slice in Fig. 3 are related to x'-y' 
coordinates of the spectral image by: 



y^y 



X X la X 



b"\ 



(2) 
(3) 



2 to a representative testing-target contained in the 
multi-PA, long-slit, wide spectral range, high-resolution 
spectroscopic survey of PNe in both hemispheres car- 
ried out with ESQ NTT+EMMI (spectral range AA 
3900-8000 A; A/AA=60 000) and Telescopio Nazionale 
Gahleo (TNG)+SARG (spectral range AA 4600-7900 A; 
A/AA=114000). 

First, a basic observative note: our long-slit spectra 
cover the entire instrumental spectral range (80 and 54 
orders for ESQ NTT-fEMMI and TNG-|~SARG, respec- 
tively), whereas long-slit echellograms of extended PNe are 
usually secured through an interference filter isolating a 
single order, to avoid the superposition of echelle orders. 
We stress that the superposition of echelle orders is unim- 
portant for an emission-line object like a PN, the actual 
limit being the effective superposition of spectral images 
in adjacent orders. Thus, a single long-slit, wide-spectral 
range echellogram provides the complete spectral image 
in a number of emissions belonging to main ionic species, 
and including Te and Ne diagnostics (|Turatto et al. 1 2002; 
IBenetti et ari 2003: and lSabbadin et al. 1 2000. 2004, 2005, 
2006). This is illustrated in Fig. 4, showing the A A 3900- 
8000 A long-slit ESQ NNT+EMMI rich-spectrum of NGC 
3918 (PNG 294.6-t-04.7. [Acker et al. , 1992). 

The order-by-order reduction procedure for long-slit, 
wide-spectral range echellograms uses IRAF and IDL 
packages, follows the standard method for bi-dimensional 
spectra (bias, CCD-flat, distortion correction, and wave- 
length and flux calibrations), and is fully described by 
ITuratto et al~l (2002). 

We select the brigh t PN NGC 7009 (the Saturn 
Nebula, PNG 037.7-34.5. lAcker et al. 1 1992) - observed at 
twelve equally-spaced PA with ESQ NTT-fEMMI (slit- 
length=60") - as a test-case for the solution of practical 
problems in the multi-dimensional turbulence-analysis. 



where a" and b" are the slopes of straight lines a' and b'. 

The foregoing analysis stresses a fundamental difference 
between imaging and spectroscopy: the former contains the 
bi-dimensional projection of the nebula on the sky, whose 
accuracy depends only on "external" factors (i. e., seeing 
and Airy disk for ground-based and orbiting telescopes, re- 
spectively), whereas the latter provides the pixel-to- pixel 
flux and velocity in different ions, and allows the nebular 
de-projection, whose accuracy depends on both "external" 
(instrumental resolution, seeing) and "intrinsic" (physical 
conditions, ionization, expansion, and turbulence) contri- 
butions. 

The odd thing is that nebular spectra - storing much 
more physical information than direct imagery - are com- 
monly integrated along the spectrograph slit, and used to 
infer the mean, average value of flux, electron temperature, 
electron density, ionic and chemical abundances, velocity, 
turbulence, etc. (thus, losing their precious spatial infor- 
mation) . 

3. The solution to practical problems in the 
multi-dimensional turbulence-analysis 

Here we solve problems related to the quantification 
and distribution of turbulent motions in real PNe. To 
this end, we apply the procedure developed in Sect. 



3.1. 1-D and 2-D turbulence-analyses of NGC 7009 close to 
its apparent minor axis 

The Saturn Nebula is an extremely complex object, con- 
sisting of several interconnected sub-systems (including 
optically-thick caps and ansae, and high-excitation polar 
streams) with different morphology, physical conditions and 
kinematics. However, nebular regions projected along and 
close to the apparent minor axis are optically thin to the UV 
radiation of the central star, and their spatio-kinematical 
structure is relatively simple and regular (jSabbadin et al."l 
2004). 

Spectral images of main ionic species present in NGC 
7009 at PA=124° (i. e. close to its apparent minor axis) 
are shown in Fig. 5, in order (left to right) of decreasing 
ionization potential, IP. They mark the spatio-kinematical 
signature in both the main shell and the outer shell of the 
Saturn Nebula. According to Sabbadin et al. 1 (2004). 

(1) the former expands at 14xp(inain shell)=4.0(±0.3) xr" 
km s~^: 



IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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Fig. 3. A model-nebular slice and the corresponding spectral image in the coUisionally-cxcited line A 5007 A of [O III] and 
in the recombination line A 6563 A of H I for different values of the instrumental spectral resolution. Left: model-nebular 
slice selected by a spectrograph slit centered on the apparent major axis of a tilted spheroid with polar diameter=20", 
axial ratio=1.4, radial thickncss(FWHM)=2.0" and yoxp(km s~^)=3.0xr". Right: corresponding [O III] (upper) and H 
I (lower) spectral images for three values of the instrumental resolution R=A/AA. 



(2) the latter is slightly blue-shifted (by 3(±1) km s'^), 
and moves at 14xp(outer shell)=3.15(±0.3)xr" km s~^; 
and 

(3) de-projected tomographic maps in different ions are 
nearly circular rings. 

Only medium- to high-excitation emissions arc contained 
in Fig. 5, NGC 7009 being optically-thin along and close 
to the apparent minor axis, with a modest concentration of 
low-ionization species (O", 0+, S+, N+ etc.). 

We start with the cspl (common to all PA, and providing 
the 1-D turbulcnce-map), and then extend the analysis to 
the whole spectral image. 

3.1.1. The 1-D turbulence-profile from the cspl 

Overall cspl flux-profiles in NGC 7009 at PA=124° are 
shown in Fig. 6. Besides the large broadening of H I and 
He II emissions, and the faint blue-shifted (red-shifted) tail 
characterizing the fine-structure profile of He II A 4686 A 
(He I A 5876 A), Figs. 5 and 6 evidence a strong line- 
asymmetry due to the slightly blue-shifted outer shell. 

Thus, we extract the net spectral contribution in both 
shells through a multi-Gaussian analysis (SPLOT routine 
of the IRAF package), as illustrated in Fig. 7 for [O HI] A 
5007 A. Such a detailed deconvolution highlights all kine- 
matical differences between nebular sub-systems, and in 
particular, the broader velocity profile of the outer shell 
iWtotai(hliic)=21.3 km s-\, and Wcoiaiircd)=20.0 km s'^ 
in Fig. 7, to be compared with the corresponding value of 



1 1 e J 1 1 > 1 1 . J- 




Fig. 6. The logarithmic cspl flux-profile in main emissions 
of NGC 7009 at PA=124°. The observed (i. c., micorrected 
for interstellar extinction) flux is in erg s^^ cm~^ px~^. The 
X-axis is centered on the systemic velocity. 

10.0 km s~^ for both blue- and red-shifted [O III] compo- 
nents of the main shell) . 

We analyze these two nebular sub-systems separately. 

3.1.1.1. The cspl of the main shell 
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1- [Nelll] 3968 A 

2- HI 4340 A 

3- [OIII] 4363 A 

4- Hell 4686 A 

5- [ArlV] 4711 A 

6- [ArlV] 4740 A 

7- HI 4861 A 

8- [OIII] 4959 A 

9- [OIII] 5007 A 

10- Hell 5411 A 

11- [NIX] 5755 A 

12- Hel 5876 A 

13- [01] 6300 A 

14- [SIII] 6312 A 

15- [01] 6363 A 

16- [ArV] 6435 A 

17- [Nil] 6548 A 

18- Hell 6560 A 

19- HI 6563 A 

20- [Nil] 6584 A 

21- [SII] 6717 A 

22- [SII] 6731 A 

23- [ArV] 7005 A 

24- [Arlll] 7135 A 

25- [Arlll] 7751 A 




Fig. 4. The upper panel shows the main part of a AA 3900-8000 A, long-slit (30") ESO NTT+EMMI spectrum (echelle 
order separation=4" (blue) to 13" (red)) of NGC 3918 at PA=10r (apparent nebular size ~27"). The brightest emissions 
are identified. The lower part of the figure is an enlargement of the red spectral region; it enhances the richness of physical 
information stored in each echellogram covering 80 orders (to be compared with the single-order output provided by the 
usual observational procedure inserting an interference filter). 
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NGC 7009 PA=124 
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HE II [Ar IV] [O III] [Ar III] He I HI 

4686 A 4711 A 5007 A 7135 A 5876 A 6563 A 

(x45) (xlSO) (xl) (xBS) (x93) (x7) 



Fig. 5. The NTT+EMMI spectral image of NGC 7009 at PA=124° in six ionic species, in order - left-to-right - of 
decreasing IP. The slit orientation is superimposed to the left-most (He II) spectral image. The original flux is multiplied 
by the factor given in parenthesis, to make each emission graphically comparable with [O III] A5007 A. More details are 
in the text. 



The instrumental spectral resolution of our echel- 
lograms, estimated from night-sky lines and comparison 
arc spectra, amounts to A/AA=60000 (i. e., Winstr=5.0 
km s~^), and is common (to within 0.5 km s~^) to all 
ionic species. Moreover, we refer to Sect. 2 for the fine- 
structure broadening (W/s) in different recombination lines 
of hydrogen- like ions. Other relevant parameters derived (or 
adopted) for the determination of turbulence, Wturb, in the 
cspl of the main sheU of NGC 7009 at PA==124° are listed 
in Table 1, where: 

- cols. 1 and 2 identify the emission line and the parent 

ion; 

- cols. 3 and 4 give the Wtotai for blue- and red-shifted 

components of the main shell; 

- col. 5 contains th e adopted value of Te (according to 

[Hyung &: AlieF] 1995a, 1995b; and ISabbadin etldTI 
2004)^ 

- cols. 6 and 7 give the FWHM of E-SE and W-NW 

zvpc-peaks in the spectral image at PA=124°, pro- 
viding - after correction for seeing-t-guiding - the ra- 
dial thickness (Ar") of each emitting layer in the 
main shcU of NGC 70090 From Ar"(cspl)~Ar"(zvpc), 
V'cxp(niain shell)=4.0xr" km s~^, and AVi2xp(niain 
shell)=Wc,rad=4x Ar" km s^^, we have 

W,ra. - 4X [( ^^'^^^-^"^+^'^^"-""^ )--Ty^,,/■^ (4) 

* In general, the value of Ar(cspl) can be inferred from the 
Ar(zvpc) trend at different PA. In the present case of NGC 7009, 
Ar(zvpc) is nearly constant along and close to the apparent 
minor axis of the nebula (i. e., at PA=4°, 124°, 139°, 154°, and 
169°). 
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Fig. 7. Multi-Gaussian analysis of the cspl velocity-profile. 
The relative [O III] A 5007 A flux-profile observed in NGC 
7009 at PA=124° (empty circles) is deconvolved in the main 
shell (long-dashed curve) and outer shell (short-dashed 
curve) contributions, and in their sum (continuous curve). 



where - the seeing-fguiding contribution given 

by the FWHM of the central star's continuum - amounts 
to 0.84(±0.05)" for afl emissions; and 
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- cols. 8 and 9 contain the resulting value of Wturb for 
blue- and red-shifted components of the cspl profile (as 
given by Eq. (1)). They provide the 1-D overall turbu- 
lence map in the main shell of NGC 7009 through the 
T4xp(km s~^)=4.0xr" relation. 

In spite of the large number of involved parameters, the 
estimated turbulence-accuracy in Table 1 is ±2 to ±4 km 
s~l for strong to faint coUisionally-excited lines, ±3 km s""'^ 
for He I A 5876 A, and ±4 km s^^ for H I Ha and He II A 
4686 A. 



3.1.1.2. The cspl of the outer shell 

The application of the foregoing procedure to the cspl 
of the outer sheU in NGC 7009 at PA=124° gives unre- 
liable turbulence values due to the huge value (and large 
uncertainty) of Wgrad, which is the predominant broaden- 
ing component of the velocity profile. 

We consider A 5007 A of [O HI] (Figs. 5, 6, and 7). On 
the one hand, we obtain W^tota/(blue)=21.3(±0.3) km s^^ 
and Wtotai (red)=20.0(±0.3) km s~^ On the other hand, the 
[O HI] thickness Ar(outer shell) - given by the zvpc profile 
corrected for seeing-|-guiding - amounts to 6.5(±0.5)"; it 
provides VFgrad=20.5(±3.5) km s~^ (through the T4xp(km 
8""'^)= 3.15(±0.3) xr" relation valid in the outer shell). 

Thus, the mere acritical, and slavish application of Eq. 
(1) to [O HI] blue- and red-shifted components of the cspl 
in the outer shell furnishes VFfnr6(blue)=-4. 5(1^9^3^) km s~^, 

and WtMr6(red)=-8. 6(^5^9^) km s^^. Even worse results are 
obtained from other, fainter emissions in Figs. 5 and 6. 

All this stresses the leading role of the spatio-kinematical 
reconstruction, whose uncertainties mark actual limits of 
the turbulence-analysis in PNe (to be noticed: in the ab- 
sence of a spatio-kinematical recovery, the broad velocity 
profile in the outer shell of NGC 7009 could be ascribed to 
turbulent motions). 

In next sections, dedicated to the 2- and 3-D mapping of 
turbulence, we focus on the net spectral-profile of the main 
shell, ignoring the outer shell (and other sub-systems) of 
the Saturn Nebula. 



model nebula 



3.1.2. Ionic and overall 2-D turbulence-maps in the main 
shell of NGC 7009 at PA=124° 

According to Sect. 2, the analysis of the velocity profile 
in all dlpl of a spectral image provides (through Eq. (1)) 
the 2-D turbulence-map of the parent ion within the entire 
nebular slice covered by the spectrograph slit. 

The [O HI] results for the main shell of NGC 7009 at 
PA=124° are graphically shown in Fig. 8, where each dlpl is 
identified by its angular distance from the cspl (d"), in units 
of the zvpc ionic radius r(0+"'")=5.55". For practical rea- 
sons. Fig. 8 contains Wtota/ (blue) and Wtotai (red) for dlpl 
at low-to-mean latitudes of the spectral image, whereas at 
high-latitudes (a) blue- and red-shifted components tend to 
merge, (b) spatial properties of the gas prevail on dynamical 
ones, and (c) the determination of Wtotai becomes problem- 
atic. Thus, for all high-latitude dlpl (i. e., d/r(O++)>0.85), 
we are forced to measure the "overall" FWHM of the de- 
convolved spectral profile. Woven corresponding to 
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Fig. 8. Main properties of dlpl in the [O HI] main shell 
of NGC 7009 at PA=124° (sampled along the slit with a 
three-pixel step, i. e., 0.84"), superimposed to the adopted 
model-nebula mimicing the spatio-kinematics of the main 
shell in the Saturn Nebula, and characterized by different 
values of Wturb. The x-axis contains the angular distance 
from the cspl (d), in units of the zvpc ionic radius r(0''"''). 
At low-to-mean latitudes (i. e., d/r(O+"'")<0.85) W refers 
to Wtotai for blue- and red-shifted components of each dlpl 
(squares), whereas at high-latitudes it provides the FWHM 
of the whole velocity-profile, i. e., Wover (circles). 

where 14 is the radial velocity of red- and blue-shifted com- 
ponents in each dlpl. 

For comparison, Fig. 8 also shows corresponding re- 
sults for the adopted model-nebula (assuming different 
values for the turbulence). This circular ring- model 
mimics [O HI] spatio-kinematical properties of the main 
shell of NGC 7009 at PA=124° (i. e., r(0++)=5.55", 
Ar(O++)corr/r(O++)=O.270 V^^p^A.Oxi" km s"! and 
To^lO^ K), and is observed with the same instrumentation 
and seeing-l- guiding conditions as NGC 7009. 

Using the T4xp(main shell)=4.0xr" km s~^ relation and 
simple geometrical considerations, we are able to map tur- 
bulent motions within the de-projected [O HI] spectral im- 
age at PA=124° (ionic 2-D turbulence-map), and the en- 
tire procedure can be repeated for other ionic species. The 
combination of all ionic 2-D tomographic-maps provides the 
overall 2-D tomographic-map of turbulence across the main 
shell of NGC 7009 at PA=124°. This is shown in Fig. 9, 
where: 

- the Ha contribution has been omitted (for obvious rea- 

sons); 

- for He II we have used the 2-D turbulence map given by A 

6560 A (this line is fainter, but also a better turbulence 
diagnostic than A 4686 A); 



Wo. 



Vrired)-Vriblue)- 



,Wtotai{blue) + Wtotaiired) 



),(5) 



^ The circular assumption applies only to untilted spectral 
images at PA=124°. In the other observed PA oi NGC 7009, 
the model-nebula is a tilted ellipse (see Fig. 6 in iSabbadin et al~l 
2004). 

« Of course, Ar(0++)co,..=(Ar(0++)L - "^1+9))°'' ■ 
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Table 1. The cspl-turbulence in the main shell of NGC 7009 at PA=124°. 
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Fig. 9. The preliminary (i. e., under-sampled) tomographic 
turbulence-map in the main shell of NGC 7009 at PA=124°, 
obtained by assembling low-to-mean latitude 2-D turbu- 
lence maps in [O III], [Ar III], [Ar IV], He I and He II (A 
6560 A). Same orientation as Fig. 5. 



- the 2-D turbulence reconstruction stops at high-latitudes 

(since the turbulence-accuracy decreases for d/r>0.85, 
due to the merging of blue- and red-shifted compo- 
nents); and 

- pixelization is due to the discrete dlpl-sampling in each 

spectral image, and fuzziness to the partial superposi- 
tion of turbulence- values in different ionic species. 



3.2. The overall 3-D turbulence-map in the main shell of 
NGC 7009 

The application of the whole procedure described in the 
foregoing section to all observed PA of NGC 7009 extracts 
the overall 2-D turbulence-map in twelve, equally spaced 
radial slices of the nebula covered by the spectrograph slit. 
Their assembly by means of our 3-D rendering program 
( [Ragazzoni et al.~| 2001), gives the spatial distribution of 



the gas-turbulence within the (almost) complete main shell 
(overall 3-D turbulence-map). 

The preliminary (i. e., under-sampled) overall 3-D 
turbulence-map in the main shell of the Saturn Nebula is 
shown in Fig. 10 and in MOVIE N. 1 (on-line material), 
where: 

- the IVt„rf)-range is 2.0 km s^^ (faintest regions) to 8.0 
km s~^ (brightest regions); 

- the left-most (00,00) panel in Fig. 10, and the fixed im- 
age at the start of the movie represent the nebula as 
seen from the Earth (North is up, and East to the left); 
and 

- the irregular, vertical, empty-band visible at (00,90) 
identifies high-latitude dlpl (i. e., d/r>0.85) in original 
spectral images. 

Although the detailed recovery and analysis of the 1- 
, 2-, and 3-D turbulence in NGC 7009 are beyond aims 
of this section (they are deferred to a dedicated paper, in 
preparation), we wish to underline the patchy turbulence- 
distribution within the main shell, and modest turbulence- 
values of innermost nebular regions at the highest ioniza- 
tion degree (i. e. He II and Ar IV emissions). In particular, 
the latter result 

- is surprising for an extended X-ray source 

([Guerrero et al. I 2002) like NGC 7009, whose hot 
and luminous powering-star (lo g T»(K)=4.95, log 
L,/Lq =3.70; ISabbadin et al. I 2004) is ejecting 
fast matter at a high - although uncertain - rate 
(terminal velocity~2700 km s~^, and mass- loss 
in the range 3.0x10"^° to 1.0x1 0"^ Mq yr"^; 
[Cerruti-Sola fc Pe rmotto] 198 5, 19 89; IBombeck et al. I 
1986; Hutscmckcr s fc Surdej ] 1989; ITinkler fc Lamers I 
2002); 

- contrasts with hydrodynamical expectations for shocked 

nebular regions interacting with the hot-bubble 
generated by fast stellar winds CFrank fc Mellem"a~l 
1994; iMellema I 19 97; IDwarkada s fc Balick] 1998; 
Garcia-Segura et al. | 1999); and 

- needs confirmation by the spatio-kinematical analysis of 

forbidden lines at an even higher ionization degree (e. 
g., A 7005 A of [Ar V], IP=59.8 eV, and A 3425 A of 
[Ne V], IP=97.1 eV ). 

4. Discussion 

Current theoretical models and radiation-hydrodynamical 
simulations of PNe predict that random, unsteady, non- 
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Fig. 10. Snapshots of the turbulence-structure in the main shell of NGC 7009 for a rotation around the North-South 
axis centered on the exciting star. The VFt^rb-range is 2.0 km s^^ (faintest regions) to 8.0 km s^^ (brightest regions). 
The line of sight is given by {9, ip)y where 9 is the zenith angle, and ip the azimuthal angle. The left-most (00,00) image 
corresponds to the nebula as seen from the Earth (North is up and East to the left). The complete 3-D turbulence- movie 
is shown in the electronic version of the paper (on-line data) . 



thermal motions at a local scale constitute a fundamental 
characteristic of the expanding gas. 

Actually, our present knowledge on the turbulence- 
structure of real PNe is nebulous, since (1) this elu- 
sive parameter has modest kinematical effects on the 
complex velocity profile, and (2) practical determina- 
tions - providing a unique, average turbulence- value 
for each PN (jGuerrero et~an 1998; lAcker et al.~l 2002; 
[Gesicki et al ~| 2003 , 2006; [Gesicki fc Zijlstra | 2003, 2007; 
and iMcdina et al. 1 2006) - are affected by one or more weak- 
enesses, namely: 

(i) poor observational material (e. g., insufhcient spectral 
resolution, very partial coverage of the nebula, limited 
number of emissions); and/or 

(ii) coarse data reduction (integration along the slit, as- 
sumption of spherical symmetry); and/or 

(iii) inadequate analysis (omission of some important 
broadening factor, absence of a spatio-kinematical re- 
construction) . 

For these reasons, we have carried out a detailed 
analysis based on the three-step methodology (spatio- 
kinematics, tomography and 3-D recovery) developed at the 
Astronomical Observatory of Padua. The disentanglement 
of all broadening components in the velocity profile (instru- 
mental resolution, thermal motions, turbulence, expansion 
gradient, and fine-structure of hydrogen- like ions) allows us 
the extraction of disordered non-thermal motions, and their 
mapping in 1-, 2-, and 3-dimensions. 

Although the observed velocity-structure of a PN 
depends on several factors - since Wgrad—i{^omc dis- 
tribution + expansion law), Wi„s(r=f (instrumental set- 
up), iyt;ier=f (electron temperature), VFturb=f (turbulence), 
and W^s=f (atomic properties) -, we argue that the final 
turbulence-accuracy mainly depends on the reliability of 
its spatio-kincmatical reconstruction. 

Thus, the detailed study of turbulence, of its stratifica- 
tion within each target and (possible) systematic variation 
among different sub-classes of PNe (A) is strictly connected 
to the general problem of de-projecting the bi- dimensional, 
apparent morphology of a three-dimensional mass of gas; 
and (B) needs (and deserves) long-slit, high-spectral res- 
olution, multi-PA echellograms over a wide spectral range 
containing emissions at low-, mean-, and high-ionization, 
to be reduced and analyzed with a straightforward and versa- 



tive methodology extracting the whole physical information 
stored in each frame at best. 

We underline that the foregoing paradigm - implying 
the recovery and pixel-to-pixel analysis of a large sample 
of emission lines in an adequate number of PA - should 
be extended to all observational parameters (e. g., flux, 
electron temperature, electron density, ionic and chemical 
abundances, etc.), due to the spatio-kinematical florilegium 
exhibite d by PNe (jTuratto et al . '2002:' Benetti et al."l 2003: 
■Sabba din et al. 1 2000. 2004, 2005, 2006). 

We end pointing out that the methodology developed 
here for PNe actually applies to all types of expanding neb- 
ulae (including nova and supernova remnants, shells around 
Population I Wolf-Rayet stars, nebulae ejected by symbi- 
otic stars, bubbles surrounding early spectral-type main se- 
quence stars, etc.) covered at adequate "relative spatial" 
and "relative spectral" resolutions. 



5. Conclusions 

So far, the quantification of stratified turbulent motions in 
the ionized gas of PNe - as predicted by hydrodynamical 
simulations and theoretical models - escaped a direct ver- 
ification because of its strict connection with the general, 
long-standing astronomical problem of de-projecting the bi- 
dimensional, apparent morphology of a three-dimensional 
mass of gas. The solution is now approached by spatio- 
kinematical, tomographic, and 3-D rendering analyses ap- 
plied to long-slit, multi-PA, high-resolution echellograms 
over a wide spectral range. 

We develop an accurate procedure disentangling all 
broadening components (instrumental resolution, thermal 
motions, turbulence, expansion gradient, and fine structure 
of hydrogen-like ions) of the velocity profile. This overcomes 
main weaknesses affecting previous determinations, and al- 
lows us to quantify and map in 1-, 2-, and 3-dimensions the 
gas-turbulence in PNe. 

The multi-dimensional turbulence-analysis - soon ap- 
plied to a representative sample of targets in different evo- 
lutionary phases, and combined with expectations from up- 
dated radiation-hydrodynamical simulations - will provide 
new, deep insights on different shaping agents of PNe (ion- 
ization, wind interaction, magnetic fields, binarity of the 
central star, etc.), and on their mutual interaction. 

Last, we emphasize all advantages - and encourage the 
adoption - of long-slit echellograms covering the entire in- 
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strumental spectral range, since the usual insertion of an 'NGC 7009' on paged 
interference filter isolating a single echelle order proves su- 
perfluous and strongly limitative in most cases... it is like 
driving a Ferrari with the hand-brake applied. 
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